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Recent magnetization and susceptibility measurements on Cu3Cl6(H20)2-2H8C4S02 by Ishii et.al. 
[J. Phys. Soc. Jpn. 69, 340 (2000)] have demonstrated the existence of a spin gap. In order to explain 
the opening of a spin gap in this copper-trimer system, Ishii et.al. have proposed a frustrated trimer 
chain model. Since the exchange constants for this model have not yet been determined, we develop a 
twelfth-order high-temperature series for the magnetic susceptibility and fit it to the experimentally 
measured one. We find that some of the coupling constants are likely to be ferromagnetic. The 
combination of several arguments does not provide any evidence for a spin gap in the parameter 
region with ferromagnetic coupling constants, but further results e.g. for the magnetization process 
are in qualitative agreement with the experimental observations. 



PACS numbers: 75.50.Ee, 75.40. Mg, 75.45. -hj 



I. INTRODUCTION 

The trimerized S = 1/2 Heisenberg chain in a strong 
external magnetic field has already received a substan- 
tial amount of theoretical attention, one reason being a 
plateau at one third of the saturation magnetization in 
the magnetization curve Some frustrated variants 

of the trimer model have also been investigated P-[lC|] 
since they can be shown to have dimer groundstates and 
thus a spin gap. 

While many materials with trimer constituents exist 
(see e.g. [0), the behavior in high magnetic fields has 
been investigated only in a few of them, for instance 
in 3CuCl2-2dioxane Also Cu3Cl6(H20)2-2H8C4S02 
belongs to the known trimer materials [y_3||lj] , but its be- 
havior in a strong magnetic field has been measured only 
recently and at the same time its magnetic suscep- 
tibility has been remeasured. Surprisingly, a spin gap of 
about 3.9Tesla (that is roughly 5.5K) is observed both in 
the magnetic susceptibility of Cu3Cl6(H20)2-2H8C4S02 
as well as in the magnetization as a function of ex- 
ternal magnetic field. This system probably exhibits also 
a plateau at one third of the saturation magnetization in 
addition to the spin gap. 




FIG. 1. The frustrated trimer chain model. All corners and 
intersections carry a spin 1 /2 coupled with exchange constants 
indicated by the connecting lines. 

Motivated by the crystal structure [|3| , the authors of 
1 15 1 have proposed the following model (see also Fig. |^) 



L/3 



H — Ji {Ssi • Ssi+i + Ssj+i • 83^+2} 

1=1 

L/3 

+ J2 ^ 83^+2 • S3i_|_3 
i=l 
L/3 

+ •^3 ^ {831+1 • 83^+3 + 83^+2 • 83^+4} 
1=1 
L 

-hY^st. (1.1) 



i=l 



Since the spin is localized on Cu^"*" ions, the 8i are spin- 
1/2 operators at site i. 



In (1.1 



the reduced field h is 
related to the physical field H by h — gfJ-sH in units 
where A;b = 1. The numerical prefactor is determined 
by Ms ~ 0.67171K/Tesla as well as the value of g which 
for the present material is slightly above 2 (the precise 
numerical value depends on the direction of the external 
magnetic field relative to the crystal axes). 

For a study of the phase diagram of the Hamiltonian 



(1.1), it is useful to observe that the Hamiltonian and 
therefore also the phase diagram are invariant under the 
exchange of Ji and J3 such that one can concentrate e.g. 
on IJ3I < |Ji|. In fact, the h — phase diagram with 
antiferromagnetic exchange constants [Ji > 0) has been 
explored in |l8[ u sing bosonization and exact diagonaliza- 
tion (see also ]19|) determining in particular a parameter 
region with a spin gap. Very recently, this was comple- 
mented by a computation of the magnetization curve at 
some values of the parameters using DMRG ||2^. The 
investigations of p"8|-p0t conc entrated on the region with 
all coupling constants in ( |l.l| ) antiferromagnetic (Ji > 0) 
because [15| suggested that this should be appropriate 
for Cu3Cl6(H20)2-2H8C4S02. However, the parameters 
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relevant to the experimental system have not really been 
determined so far. We believe that this is an important 
issue in particular in view of the fact that according to 
the crystallographic data jl^ , all angles of the Cu-Cl- 
Cu bonds lie in the region of 91° to 96° - a region where 
usually no safe inference on the coupling constants Ji can 
be made, not even about their signs. We will therefore 
develop a high-temperature series for the magnetic sus- 
ceptibility of the model (1.1) and use it to determine the 
coupling constants from the experimental data |p^ . It 
will turn out that some coupling constants are likely to 
be ferromagnetic^ i.e. the experimentally relevant cou- 
pling constants lie presumably outside the region studied 
so far. We then proceed to study more general properties 
of the model and to address the question of a spin gap in 
the relevant parameter region. We use mainly perturba- 
tive arguments supplemented by numerical methods. 

Some supplementary results on the trimer model are 
contained in the appendices or can be found in pl|. 



II. MAGNETIC SUSCEPTIBILITY AND 
SPECIFIC HEAT 

A. High-temperature series for zero field 

First we discuss some high-temperature series in zero 
magnetic field. We have used an elementary approach to 



perform the computations. Denote the Hamiltonian of a 
length L chain with h = Qhy Hq. Then the fundamental 
ingredient for any higher-order expansion is that con- 
tributions of tr {H^^ to suitable physical quantities be- 
come independent of the system size L if one uses a long 
enough chain with periodic boundary conditions. The 
concrete Hamiltonian Hq given by (1.1) must be applied 
2L/3 times to wind once around the system and to feel 
that it is finite. On the other hand, contributions from 
tr {H^^ with N < 2L/3 are independent of L. We have 
used this observation to determine the high-temperature 
series by simply computing the traces for the lowest pow- 
ers iV on a chain with a fixed L and periodic boundary 
conditions j2^. Just two small refinements to this ele- 
mentary approach have been made. The first one is that 
we computed the traces separately for all subspaces of 
the z-component of the total spin S^^^. This is already 
sufficient to obtain series for the specific heat Cy and the 
magnetic susceptibility x- The second one is to make also 
the order 2L/3 usable: At this order, only the coefficient 

L /3 L /3 

of Jj^ ' Jg ' is affected by the finiteness of the chain and 
this coefficient can be corrected by hand using results for 
a Heisenberg ring of length 2L/3. 

For notational convenience, we introduce the partition 
function for L sites by 



Zr 



tr (e 



(2.1) 



with ksT = 1/(3. 



The lowest orders of a reduced magnetic susceptibility x are found to be 



Xrod.(/3) 



1 



(3LZl dh^ 
4 24 



tr (e-'^t^^ 



h=0 



tr 



((^tot) 



2 g-/3i/o 



Zr 



0' 



(2J3 + J2 + 2Ji) - ^ (J2 + 4 + Jl - 6J3J1 - 2J2J1 - 2J3 J2) 



96 



1152 
4608 



(8 J| - 3Jp2 + 6J| Jl + 6J3 J| + 8Jf - 3J2 Jl + Jl) 
(-28J3V2 Jl - 28J3 Jl J2 + 36J| Jf + 8Jpi + 8Jp^ - 2J| Jl - 34J3 Jf - 34J| J; 



-IOJ2J1 - 2J3 J| - 10J|J2 - 28J3 Jl J| + 14Ji^ + 14J| + 5J|) 

Similarly, we obtain the lowest orders of the high-temperature series for the specific heat 

UP) P^ 



(2.2) 



16 ^ ^1 



ln(Zi 



2^1) + § (2^1 



2Jl 



- 6J3 Jl J2 



256 

+0 {13'' 



(8J3V2 Jl -t- 8J3 Jl J2 -f UJpl + 8Jp^ + 8Jpl + &jf + 6JI + J| -t- 8J3 Jl J|) 



(2.3) 



Complete 12th order versions of both series can be accessed via |2l|| . 

For a uniform Heisenberg chain (Ji = J2 and J3 = 0), the coefficients of the series for x a-nd In (Zl) /L (or Cy) 
agree with those given for instance in 123] when they overlap. 
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B. Fit to the experimental susceptibility 



Now we use our 12th order series for the susceptibility 
( ^.2[ ) to fit the experimental data and thus extract 
values for the coupling constants Ji. We used the data 
for the single crystal {H \\ 6-axis) and the polycrystalline 
sample |15| as well as some unpublished new measure- 
ments for all three axes of a single crystal |Q . For the 
polycrystalline case the average g-factor is known to be 
(7av ~ 2.1 from ESR while in the single crystal case p^ ] 
we used the g-factor as a fitting parameter. The follow- 



ing prefactors are used to match the series (2^) to the 
experimental data: 



Xoxp. {T) = 



Xrcd. [ rp 



(2.4) 



We performed fits in various intervals of temperature 
with a lower boundary (T/) lying between 150K and 
250K, while the upper boundary was kept fixed at 300K. 
Fits were performed with the raw 12th order series. For 
both experimental data sets of jl^ we obtained rea- 
sonable, though volatile fits around T/ — 150K — 250K 
yielding the following estimates: Ji = — 250K ± 40K, 
J2 = 250K±40K, J3 = -40K±30K. For the single crys- 
tal sample we additionally determined gb = 1.95 ± 0.05. 

We have further performed fits to unpublished single- 
crystal data sets where the ^-factors ar e kn own from 
ESR (2|]. When a constant is added to dJ), the data 
for all three crystal axes can be fitted consistently with 
Ji « -300K, J2 « 280K and J3 « -60K in an interval of 
high temperatures (220K < T < 300K). This set of cou- 
pling constants is in agreement with our earlier fits and 
we will use the latter in the further discussion below. 

Fig. |2| shows the measured susceptibility for the poly- 
crystalline sample jl^] together with the series result. 
Since the parameters were obtained from a fit which 
was performed with a different data set, we have used 
g = 2.03 (which differs slightly from the experimentally 
found ^av ~ 2.1) in order to obtain agreement of the raw 
series with the experimental data for T > 240K. Clearly, 
the raw series should not be trusted down into the region 
of the maximum of x where Fade approximants should be 
used instead. The region below the maximum cannot be 
expected to be described with a high-temperature series. 
The overall agreement is reasonable though the theoret- 
ical result reproduces the experimental one in the vicin- 
ity of the maximum only qualitatively. This discrepancy 
might be due to the frustration in the model which leads 
to cancellations in the coefficients. Note also that, due to 
the frustration, the maximum of x is located at a lower 
temperature than would be expected for a non-frustrated 
model with coupling constants of the same order of mag- 
nitude. Consequently, higher orders are important in the 
entire temperature range covered by Fig. ^, precluding in 
particular the analysis of the high-temperature tail of x 
in terms of a simple Curie- Weiss law. 



The agreement for intermediate temperatures can be 
improved if the maximum is included in the fitting re- 
gion and Fade approximants are used in the fit. The 
main change with respect to the fits discussed above is 
that J3 tends to be closer to Ji. However, it will become 
clear from the discussion in later sections that the region 
with J3 close to Ji is not appropriate to describe the ex- 
perimental observations of the low-temperature region. 
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FIG. 2. Experimental results for the susceptibility ('-I-') 
in comparison with the fit Ji = — 300K, J2 = 280K and 
J3 = — 60K. We show the raw 12th order series (dotted line) 
as well as several Fade approximants: [7,6] (full fine), [6,6] 
(long dashes) and [6,5] (short dashes). 

Although we are not able to determine the coupling 
constants to high accuracy, all our fits lead to the con- 
clusion that J2 should be antiferromagnetic and Ji and 
J3 (or at least one of them) must be ferromagnetic if 
one wants to model the susceptibility measured at hig h 
temperatures [T5| ] with the frustrated trimer chain (1.1). 
In view of earlier theoretical investigations [p^-pO[, this 
conclusion is somewhat surprising. Note that none of 
our fits converged to all J; > 0. Additional assumptions 
(including a constraint on the Ji) are necessary to deter- 
mine from x(T) what the optimal values of the Ji would 
be in this antiferromagnetic region and thus allow for a 
comparison with Fig. Such a fit and a comparison with 
the present one is discussed in appendix H. The upshot 
is that the experimentally observed x(^) lllSt cannot be 
explained with only antiferromagnetic Ji. 

The findings of this section necessitate a detailed re- 
analysis of the Hamiltonian (^]|) since earlier works did 
not look at the appropriate parameter region. 



III. LANCZOS RESULTS 

In order to study the zero-temperature behavior of the 
frustrated trimer chain we have performed Lanczos di- 
agonalizations of small clusters with periodic boundary 
conditions. Although computations were performed for 
various values of the parameters, we will present explicit 
results only for the final parameter set determined above. 
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Further results in the region > are in agreement with 
and are used in appendix 0. 
Fig. ^ presents the zero-temperature magnetization 
curve for the trimer chain model. Here and below the 
magnetization (M) is normalized to saturation values ±1. 
First, it is reassuring that the system still has antifer- 
romagnetic features despite two ferromagnetic coupling 
constants (note that we are now probing a region far from 
that used for determining the Ji). Since experiments 
found a spin gap [ p^ , an important question clearly is if 
we also obtain a gap from the model with these parame- 
ters. We have therefore performed a finite-size analysis of 
the gap to 5'^ = 1 excitations (corresponding to the first 
step in the finite-size magnetization curves of Fig. |^) . All 
our approaches led to results compatible with a vanish- 
ing gap. However, it is difficult to reliably exclude a gap 
of a few K with system sizes L < 30. We will therefore 
return to this issue later and assume for the extrapolated 
thick line in Fig. ^ a vanishing spin gap. In general, this 
extrapolation was obtained by connecting the mid-points 
of the steps of the L — 30 magnetization curve, except 
for (M) = 1 and (A/) = 1/3 where the corners were used. 
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FIG. 3. Magnetization curve for Ji = -300K, J2 = 280K 
and J3 — — 60K. The thick line is an extrapolation whereas 
thin lines are for finite system sizes: L = 12 (dotted), L — 18 
(short dashes), L — 24 (long dashes) and L = 30 (full). 

For the parameters of Fig. ||, (M) = 1/3 is reached 
with a magnetic field H = 20 — 25Tesla. The order of 
magnitude agrees with the experimental finding p5| ] even 
if the value found within the model is a factor of two to 
three below the experimental one. Above this field. Fig. 
^ exhibits a clear (AI) = 1/3 plateau which is expected 
on general grounds [Q^. 

We conclude this section by presenting in Fig. ^ the 
lowest three excitations for the 5^ = 1 sector as a func- 
tion of momentum k, where k is measured with respect 
to the groundstate, i.e. k — ks'=i — kos- This spectrum 
is very similar to that of an S" = 1 /2 Heisenberg chain of 
length L/3 with coupling constant J^s. ~ 16K. In partic- 
ular, one can recognize the two-spinon scattering contin- 
uum and a few higher excitations. This identification of 
the low-energy excitations of the frustrated trimer chain 
with an effective 5=1/2 Heisenberg chain is one of the 



numerical indications for the absence of a spin gap. 
60 




FIG. 4. Lowest three excitations in the — 1 sector for 
L = 24 ('+') and L = 30 ('x') with Ji = -300K, J2 = 280K 
and J3 = -60K. 



IV. THE LINE Ji = J3 



The Lanczos results of the previous section raise the 
question if the trimer chain model has a spin gap in the 
region with Ji, J3 < 0: In this region, the model behaves 
like an antiferromagnet (see e.g. Fig. ^) which is frus- 
trated since the number of antiferromagnetic coupling 
constants around a triangle is odd. Therefore, a spin gap 
appears possible in principle and we proceed with further 
arguments to decide whether it appears in the relevant 
parameter region. 

Evidence for a spin gap in the parameter region Ji, 
J3 > was actually first obtained on the line Ji = J3 
]lO| . The reason is presumably that the line Ji = J3 
can be treated analytically at least to some extent be- 
cause then the total spin is locally conserved on each 
bond coupled by J2. In fact, one can easily discuss the 
entire magnetization process |^5| and not just the ques- 
tion of a spin gap and we refer the interested reader to 
^l| for some comments on this aspect. 

Recall that the model (^]^) with h = gives rise 
to three types of groundstates in different regions with 
Ji — J3 > (we will assume J2 > throughout this sec- 
tion) |l^: That of the S = 1/2- S = 1 ferrimagnetic chain 
for J2 < O.9O8I6J1, a spontaneously dimerized state for 
O.9O8I6J1 < J2 < 2Ji and, finally, singlets are formed 
on all bonds coupled by J2 for J2 > 2 Ji with effectively 
free spins in between. In the first and the third case, one 
finds ferrimagnetic behavior with a spontaneous magne- 
tization (M) = 1/3 and only the second region exhibits 
the requested gap. 

For Ji = J3 < 0, we found only two regions: 



1. For 



J2 > -Ji 



(4.1) 
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the groundstate is formed by singlets on the J2- 
bonds and free S = \/2 spins in between. This 
again gives rise to ferrimagnetic behavior with a 
spontaneous magnetization {M) = 1/3. 



2. When 



J2 < -Ji 



(4.2) 



the entire system behaves hke a ferromagnet. In 
this case the system is spontaneously completely 
polarized {{M) = 1). 

We conclude that -unlike for = J3 > 0- the 
groundstate is always gapless for = J3 < 
which we have argued to be more appropriate for 

CU3C16(H20)2-2H8C4S02. 



V. EFFECTIVE HAMILTONIANS FOR THE 
GROUNDSTATE 

Th e low-energy behavior of the model Hamiltonian 
(1.1) can be analyzed further using degenerate perturba- 
tion theory. Truncation at a certain order of the coupling 
constants leads to effective Hamiltonians which in some 
cases turn out to be well-known models. 

We will use the abbreviations 



"^^ - Ji ' 



- J2 ■ 



(5.1) 



A. Ji large and antiferromagnetic 

To test the method, we first consider the case of anti- 
ferromagnetic Ji. For this purpose we extend the first- 
order effective Hamiltonian of j20| for the case Ji ^ J2, 
J3 > to second order. For Ji large and antiferromag- 
netic, the groundstate-space of a trimer is given by an 
S = 1/2 representation. In this subspace of doublets, 
the effective Hamiltonian has the form of a J1-J2 chain 
when truncated after the second order: 



Ji Si ■ Si+i + J2 Si ■ s 



i+2 ■ 



(5.2) 



Here, the Si are effective spin-1/2 operators. The effec- 
tive exchange constants are 



Ji 



(-./3 + J2) 



and 



J2 
Ji 



91 ~, 



— J2 

405 ^ 



22 ~ 



135 1620 ^ 



10 ~, 



486 



243 



243 



(5.3) 



(5.4) 



In this approximation, the ferrimagnetic phase found in 
ITsf is given by an effective ferromagnetic Hamiltonian 
{Ji < 0) while the antiferromagnetic phase corresponds 
to iTi > 0. The transition line can thus be determined 
from Ji = pi. We find 



~ _ 12 ~ 
= J2 



90 



1 



— + a/17245J,2 

79 158 V 2 



3^^2^+0(i2^), 



48240 J2 + 32400 



(5.5) 



which improves the agreement of the approximation J3 = 
J2 with the numerical results of 18|. 

The dimer phase with a spin gap is characterized by 
J2/J1 > 0.241167(5) (see ||2^ and references therein). 
Using ( |5.3| ) and (^), it is found to open at J2 ~ 3.60, 
J3 « 1.361 with a square-root like behavior of J3 as a 
function of J2. Since this is not in the weak-coupling re- 
gion, it is not surprising that the numbers differ substan- 
tially from those obtained numerically in jl^ . However, 
the topology of the groundstate phase diagram comes out 
correctly from our effective Hamiltonian: In particular, 
the dimerized spin gap phase is located inside the anti- 
ferromagnetic phase and arises because of a sufficiently 
large effective second neighbor frustration J/2- 



B. J2 large and antiferromagnetic 

The preceding argumentation is not applicable to the 
region J2 > 0, Ji, J3 < 0. However, a similar case has 
been discussed earlier and J2 ^ | Ji | , | J3 1 has been 
found to be a useful limiting case. We will now analyse 
this region in the same manner as above. 

For J2 » I Ji I , I J3 1 , the spins on all J2-bonds couple 
to singlets and only the intermediate spins contribute to 
the low-energy excitations. In the space of these inter- 
mediate spins, we ca n ag ain map the Hamiltonian (1.1) 
to the Hamiltonian (5.2) to the lowest orders in Ji, J3. 
Up to fifth order, we find the effective coupling constants 
to be given by |2q] 



Ji 
J2 



1 3 ( Ji + J3 



3Ji J3 



12 4 

(Ji + J3) (107 {.11 + J|) - 4O6J1 J3 
64 



and 



J2 



(5.6) 



(5.7) 



This mapping is now applicable regardless of the sign of 
Ji and J3 as long as J2 > 0. First we consider the case 
of antiferromagnetic Ji, J3 > 0. Then the effective cou- 
pling constants are essentially always antiferromagnetic, 
i.e. Ji,J2 > leading to a frustrated chain. If Ji and 



5 



Ja are large enough, J2/J1 can exceed the critical value 
of about 0.241 (see above) and a spin gap opens. These 
observations are again in qualitative agreement with the 
phase diagram of [^sj. As for the preceding limit, one 
should not expect good quantitative agreement since the 
required values of Ji and J3 are not small but of the same 
order as J2. 

Now we turn to the mor e int eresting case Ji, J3 < 0. 
Then the oupling constant (5/7) is always ferromagnetic: 
J2 < 0. If |Ji| and IJ3I are large enough, Ji also be- 
comes ferromagnetic. This is compatible with the behav- 
ior found in section IV on the line Ji = J3 < 0. If | Ji| and 
I J3I are small, J7i remains antiferromagnetic. Since J2 is 
always ferromagnetic, no frustration arises in the effec- 
tive model and a spin gap is not expected to open. This 
is true to the order which we have considered. Higher 
orders might actually yield frustrating contributions. In 
any case, frustration is substantially weaker for ferromag- 
netic Ji, J3 < than for antiferromagnetic Ji, J3 > 0. It 
is therefore plausible that a spin gap is absent in the fer- 
romagnetic region (unless | Ji| and/or | J3I are very large 
and the present argument is not applicable). 

It should be noted that ( |5.6| ) and (^) turn out to be 
small if Ji — J3 is small. In fact, one can argue that 
the results of this section remain qualitatively correct 
for Ji — J3 small even if Ji and J3 are not separately 
small: For Ji = J3, the intermediate spins are effec- 
tively decoupled due to the presence of the singlets on 
the J2-bonds (see section |v|) . A small detuning Ji ^ J3 
generates an effective coupling of the intermediate spins 
via higher-order processes. However, the effective cou- 
pling will stay small as long as Ji — J3 is small. If one 
wants to model Cu3Cl6(H20)2-2H8C4S02, \ Ji — J-i\ must 
therefore at least be on the same scale as e.g. the field 
h « 80K required to polarize the intermediate spins lead- 
ing to (M) = 1/3 Q. This observation rules out a Ji 
very close to J3. 

Finally, we also calculated the effective Hamiltonian for 
a strong ferromagnetic intra- trimer interaction Ji. The 
problem then maps to a frustrated S' = 3/2 chain with 
four-spin interactions. Even if this is not a well-known 
Hamiltonian and the issue of a spin gap thus remains un- 
clear in this case, we present it in appendix^ in order to 
open the way for further investigation of this limit. 



VI. MAGNETIZATION PLATEAUX 

We complete our theoretical analysis with a discussion 
of plateaux in the magnetization curves of the frustrated 
trimer chain model. 

A plateau with (M) = 1/3 is abundant in the magne- 
tization curve (compare Fig. H) and can be easily under- 
stood in the limits | J2I, | J3I ^ Ji or |Ji|, IJ3I <C J2. This 
is readily done by adding the coupling J3 to the series 
of [Q. More details as well as the explicit series for the 
boundaries of the {M) = 1/3 plateau are available under 



1 21| . Here we just mention that the main conclusions of 
I p^regarding this plateau remain qualitatively unchanged 
in the presence of the additional coupling J3. 

Regarding plateaux with (M) 7^ 1/3, observe first that, 
when a spin gap opens in the frustrated trimer model, 
the groundstate is dimerized, i.e. translational invariance 
is spontaneously broken by a period two. Spontaneous 
breaking of translational invariance by a period two also 
permits the appearance of a plateau with {M) = 2/3 (see 
||29| and references therein) . We will now investigate this 
possibility further. 

First we consider the case Ji > and start in the limit 
of strong trimerization ( J2 = 0, J3 = 0). When one ap- 
plies a magnetic field hc — ^Ji, the two states |tTT) £^nd 
^(liTT) -2 ITiT) +lTTi)) are degenerate in energy. This 
degeneracy is then lifted by the couplings J2 , <^3 • The ef- 
fective Hamiltonian to first order is an XXZ chain in a 
magnetic field ||30|- |35| , p^ . We obtain the following effec- 
tive couplings for the XXZ chain: 



Jxy 



■7:^2 — -;:J'i 
6 3 

1 



J. = 3g(J2+8J3) 

/loff = h- he - ^{5J2 + 22 J3) 
3o 



(6.1) 



and therefore the effective anisotropy Acff = Jz/\Jxy \ is 



Acff = 



J2 + 8 J3 
I6J2 -24J3I 



(6.2) 



For 5/32 < J3/J2 < 7/16, we have A^g > 1 and thus a 
gap, i.e. an (M) — 2/3 plateau in the original model. A 
plateau with (M) =2/3 can be indeed observed numer- 
ically somewhere in this region (see e.g. |^). The line 
J3/J2 — 1/4 describes the Ising limit A^s = 00. 

In order to address the region of ferromagnetic Ji, 
we now start from the limit Ji = J3 = and apply a 
magnetic field he — J2. Then the two states jft) and 
(lit) — ITi)) on the J2-dimer become degenerate in 
energy while the intermediate spins are already polar- 
ized. This can be again treated by degenerate perturba- 
tion theory in 1/J2- Up to third order we find an XXZ 
chain with 



Jxy 

J2 

that is 



J1 + J3) {Ji'hy 



h I 
T2-'~2^'^ 



J1 + J3 



J1 + J3 



(6.3) 



(6.4) 



In the region where this treatment is valid, we always 
have a small Aes, i.e. no plateau at (M) = 2/3. Indeed, 
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one can see that the dimer excitations can hop at second 
order in I/J2 while up to this order all diagonal terms 
involve only a single dimer site. Thus, up to second or- 
der the diagonal terms contribute only to hce and to this 
order one obtains an XY chain in a magnetic field. A 
small anisotropy is restored at third order before terms 
that are not described by a simple XX Z chain arise at 
fourth order. 



VII. CONCLUSIONS 



We have studied the frustrated trimer chain (1.1) 
(Fig. |l|) using a variety of methods. First, we have com- 
puted 12th-order high-temperature series for the suscep- 
tibility X and specific heat. Fits of the high-temperate 
tail of the susceptibility computed from the model to the 
one measured on Cu3Cl6(H20)2-2H8C4S02 [|l5| lead to 
J2 = 250K±40K and ferromagnetic Ji = -260K±50K, 
J3 = -40K ± 30K (we showed in appendix g that xiT) 
cannot be fitted with the antiferromagnetic parameters 
proposed in p^-pO[|). We assumed that these parame- 
ters remain valid down to low temperatures since we are 
not aware of any indication of a drastic change in the 
magnetic behavior of Cu3Cl6(H20)2-2H8C4S02 as tem- 
perature is lowered. In fact, features of other experi- 
mental observations at intermediate and low tempera- 
tures are roughly reproduced with the aforementioned 
parameters: We find a maximum in x(^) in the region 
50K < T < lOOK and a smooth increase of the low- 
temperature magnetization (M) from to 1/3 as the ex- 
ternal magnetic field is increased from zero to several 
ten Tesla. From a quantitative point of view, the agree- 
ment may however not yet be entirely satisfactory: Devi- 
ations between the measured susceptibility from the one 
obtained within the model can be seen in the interval 
80K <T < 200K and the model predicts an (M) = 1/3 
magnetization for a magnetic field that is a factor two to 
three below the one actually required in the experiment. 

Probably the most exciting experimental observation 
H for Cu3Cl6(H20)2-2H8C4S02 is the existence of a 
spin gap of about 5.5K. We have therefore searched for a 
spin gap in the region of ferromagnetic Ji and J3 using 
several methods. Neither Lanczos diagonalization, dis- 
cussion of the line Ji = J3 nor an effective Hamiltonian 
for large J2 provide any evidence in favor of a spin gap 
in this parameter region. A further careful analysis of 
this issue would certainly be desireable in particular in 
view of the small size of the actually observed gap. At 
present, however, it seems likely that the model does not 
reproduce a spin gap in the relevant parameter region. 

It should be noted that the coupling constants which 
we have determined are about two orders of magnitude 
larger than the experimentally observed gap. There- 
fore, a small modification of the model is sufficient to 
produce a gap of this magnitude. The possibilities in- 
clude dimerization of the coupling constants, exchange 



anisotropy as well as additional couplings. A modifica- 
tion of the model along these lines may also help to im- 
prove the quantitative agreement with the features ob- 
served in Cu3Cl6(H20)2-2Il8C4S02 at energy scales of 
about lOOK. Further measurements are however needed 
to discriminate between these possibilities. For example, 
it would be interesting to measure the specific heat and 
compare it with our series ( ^.3| ). It emerges also from our 
analysis that a temperature of 300K is still too small to 
allow for application of a simple Curie- Weiss law to the 
magnetic susceptibility x- It would therefore be useful 
to measure x to higher temperatures in order to permit 
analysis via truncation of (2.2) after the order which 
would provide a more direct check that 2(Ji + J3) + J2 
is negative. 

However, inelastic neutron scattering would presum- 
ably be most helpful: First, this should clearly de- 
cide if Cu3Cl6(H20)2-2Il8C4S02 is really quasi-one- 
dimensional and secondly it would yield direct informa- 
tion on the excitation spectrum which could hopefully 
be interpreted in terms of coupling constants. Such a de- 
termination of the coupling constants would also circum- 
vent the question whether model parameters change as a 
function of temperature since neutron scattering is car- 
ried out at low temperatures, i.e. the temperature scale 
of interest. We therefore hope that neutron scattering 
can indeed be performed and are curious if excitations 
will be observed that are similar to those computed in 
the trimer chain model (Fig. 

The frustrated trimer chain model is also interesting in 
its own right: It has a rich phase diagram which among 
others includes many aspects of the J1-J2 chain such as 
a frustration-induced spin gap in some parameter region 
p8|-p0[. Also plateaux in the magnetization curve exist 
in this model: A plateau with (M) = 1/3 is abundant 
both in the regions with antiferromagnetic and ferrimag- 
netic h = groundstates. Also a plateau with (M) = 2/3 
can be shown to exist in the region with Ji , J3 > (see 
pot and section ^) . Like in the case of the spin gap, the 
opening of the latter plateau is accompanied by sponta- 
neous breaking of translational invariance in the ground- 
state. Amusingly, however, the (M) = 2/3 plateau opens 
already for J2, J3 <JC Ji - a region where the spin gap is 
absent. In this context of magnetization plateaux, we 
hope that the magnetization measurements can be 
extended to slightly higher fields which should unveil the 
lower edge of the (M) = 1/3 plateau. 
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APPENDIX A: ANTIFERROMAGNETIC 
COUPLING CONSTANTS 

In this appendix we discuss a fit of the magnetic sus- 
ceptibility X with antiferromagnetic coupling constants 
Jj > 0. A number of assumptions are necessary in order 
to obtain at all a convergent fit with parameters in the 
antiferromagentic spin gap region jl^-g^l- 

First we fix the ratio of the magnetic field h{{M) = 
1/3) to the spin gap hc{{M) = 0) approximately to the 
experimental value Pq] 



HiM) = 1/3) 
hc{{M) = 0) 



= 14.1, 



(Al) 



To this end, we used numerical data for h{{M) = 1/3) 
and hc{{M) = 0) on systems of size L = 12, 18 and 
24. This data was extrapolated to i = oo in the same 
manner as in p9| ], i.e. with a polynomial fit hL{{M)) — 
hoo{{M)) + a/L 6/L^. For the spin gap, this amounts 
to reproducing t he computation of |]l^. The numerical 
solutions to eq. (|Al|) were then approximated by 



Js = 0.3 J, - 0.85 



0.63 



(A2) 



where we used the notation (5.1). An analytic for mula 
was needed i n o rder to im plement the constraint (Al) 
by inserting (A2) into (2.2) before performing a fit. Eq. 



(g|) is vahd for 0.6 < J2 < 2. 

The constraint ( |A2| ) is still not sufficient to ensure an- 
tiferromagnetic Ji > with 0.6 < J2 < 2. To achieve 
this goal, we had to make the foll owing further adjust- 
ments when fitting our series ( |2.2| ) to the experimental 
data |15|: 



1. Keep g as a fitting parameter. 



2. add a constant to (2^) and use this as another pa- 
rameter in the fit, 

3. start fitting at low temperatures Ti k, lOQK. 

Note that both g and the additive constant turn out to be 
quite large. For example, for the parameters used in Fig. 
|5|, we found g « 2.9 and an additive constant of about 
—0.1 7 • lO^'^K^^/Zcs. This means that the prefactor in 
( ^.4|) is off by a factor of about 2 from the value deter- 
mined by ESR and that the absolute value of the additive 
constant is almost 40% of the susceptibility observed at 
T = 300K ! 

On the basis of these unrealistic parameters, one could 
already discard this fit to x(^)- Nevertheless, we com- 
pare it to the one shown in Fig. Fig. ^ shows the mea- 
sured susceptibility for the polycrystalline sample [|l5| to- 
gether with the series evaluated at Ji — 120K, J2 = 141K 



and J3 = 79K. This parameter set is close to parame- 
ters proposed in jl^. This proposal was based on two 
assumptions: 1) The model should give rise to the ex- 
perimentally observed spin gap of around 5K jl^. 2) 
The maximum of x is located at T « 0.7Ji. While we 
do indeed reproduce the spin gap rather accurately, the 
second assumption is falsified by our computation: The 
frustration pushes the maximum of x(^) aga-in to lower 
temperatures as compared to a non-frustrated system. 

1 



o 



0.75 



0.5 



0.25 




i 





50 100 150 200 250 300 
T[K] 

FIG. 5. Experimental results for the susceptibility ('+') in 
comparison with the fit Ji — 120K, J2 = 141K and J3 — 79K. 
We show the raw 12th order series (dotted line) as well as the 
[7,6] Pade approximant (full line). 

Fig. H should be compared to Fig. ^. The seemingly 
better agreement in the region lOOK < T < 200K is due 
to including this temperature interval in the fit for Fig. 
^, but not in Fig. ^j. Note that the \ Ji\ are now smaller 
by a factor of about two than those used in Fig. ^. One 
would therefore expect better convergence in the vicinity 
of the maximum of x{T), i.e. for 50K < T < lOOK. This 
expectation is confirmed by the fact that in Fig. ||, the 
[7,6], [6,6] and [6,5] Pade approximants are indistinguish- 
able. However, while the series reproduces the maximum 
roughly in Fig. ^ this is definitely not the case in Fig. ||. 
The better agreement of the fit in Fig. || with the experi- 
mental data at r « 70K is particularly remarkable since 
this temperature range is far from the fitting region in 
this case, while closeby temperatures were used in Fig. |^. 
In combination with the unrealistic assumptions needed 
to obtain a convergent fit with all Ji > one can there- 
fore conclude safely that only antiferromagnet coupling 
constants are not suitable for describing the experimental 
data |15| for the susceptibility x(T). 



APPENDIX B: EFFECTIVE HAMILTONIAN FOR 
Ji LARGE AND FERROMAGNETIC 

For a strong ferromagnetic intra-trimer interaction Ji , 
the noninteracting groundstates are built from products 
of trimer 5 = 3/2 states. Up to second order we find 
the following effective Hamiltonian in this subspace of 
low- lying trimer quartets: 
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^cff. — Ja ■ S^+i + J\y Si ■ Si+2 

i i 

+Jc5I(s«-s.+i)' 



-4- ' Sj+i)(Si+i ■ s 



i+2j 



+ (Si+2 • Si+i)(Si+i • Si)} , 

where the are now effective spin-3/2 operators. 
The couphng constants are found to be: 



(Bl) 



Ja-^(J2 + 2J3) + 



197J| + 212 J2 J3 + 212 J| 



2592 1 Ji I 



Jc 

Jd 



2Jl + 5J2 J3 + 2J| 

41J| + IOOJ2 J3 + 36J| 
1296 1 Ji I 
4(2 J| + 5 J2 J3 + 2 J|) 
243|Ji| 



(B2) 



Even if this effective Hamiltonian is not a well-known one, 
it is clear that there is no spin gap in first order, since 
then the system is effectively a nearest neighbor S — i/2 
Heisenberg chain which is either gapless ( J2 + 2 J3 > 0) 
or ferromagnetic (J2 + 2J3 < 0). 

If one neglects the Jc and Jd terms, one obtains a frus- 
trated S = 3/2 chain which has been investigated with 
DMRG and leads to a gap for Jb/Ja > 0.3 @. It seems 
to be possible to obtain antiferromagnetic Ja and Jb in 
this region if J2 and J3 are chosen suitably and large 
(a regio n inc luding the coupling constants determined in 
section II B ). However, then one is not in the pertur- 
bative region anymore and the Jc and Jd terms may 
also become important. Further discussion is therefore 
needed for reliable conclusions about a gap on th e ba sis 
of the Hamiltonian (Bl) with coupling constants 
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